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SUMMARY 


The internal conversion process has been analyzed by taking advantage of the simplifications 
in the mathematical formalism which are possible near the conversion threshold. The results 
are applied to M1 conversion in the K-shell. Corrections for various nuclear charge distri- 
butions can be easily made and a general insight into the nature of the problem can be 
obtained, allowing such predictions as the maximum deviation possible from the so-called 
normal values for the internal conversion coefficients that were first computed by Sliv and 
Band. 


Introduction 


In nuclear spectroscopy, the study of the internal conversion of nuclear 
gamma-rays supplies important information about the excited states of nuclei. Dur- 
ing the last fifteen years there has been a rapid improvement in the facilities 
for making more precise calculations and accurate experimental determinations 
of the internal conversion coefficients. The development of fast electronic com- 
puters has made their accurate computation possible (Rose et al, 1951, Sliv and 
Band 1956 and 1958, Rose 1958) and the beta-spectrometrie method for the ex- 
perimental study of the internal conversion in complex spectra has been im- 
proved (Hultberg and Stockendal 1958). 

With the publication of extensive tables for the K- and L-shells (Sliv and 
Band 1956 and 1958, Rose 1958) and the M-shell (Rose 1958), the experimen- 
talist has obtained a fairly easy access to the internal conversion coefficients 
for such values of the gamma-ray energy and the multipolarity which are of 
interest in practice, at least for the K- and L-shells. 

In the earlier calculations of Rose, the nuclear dimension were neglected. 
This is generally a very good approximation but, in spite of the small dimen- 
sions of the nucleus as compared with the atomic dimensions and the wave 
length of the nuclear gamma-ray, the nucleus can have a considerable influence 
on the internal conversion rate. This circumstance was first realized by Sliv. 
(Sliv 1951, Sliv and Listengarten 1952.) In 1956, Sliv and Band published com- 
prehensive tables of conversion coefficients. These computations were based on 
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the assumption that the displacement of nuclear matter takes place on the 
nuclear surface only and this means that there is no interaction between nuclear 
matter and the atomic shell electrons when the latter are inside the nucleus. 
Mathematically this. corresponds to cutting off the radial part of the transition 
matrix elements inside the nucleus. The presence of the nucleus makes it nec- 
essary to correct the electron wave functions. The values for the internal con- 
version coefficients obtained using this model differ by as much as 50 per cent 
from the point-nucleus values which were computed by Rose. These large dif- 
ferences are mainly caused by the removal of the singularities in the radia- 
tion field as well as in the electron wave functions for certain states as a con- 
sequence of the finite nuclear size. As soon as the singularities are removed, 
the internal conversion coefficients are rather insensitive for details in the nuclear 
structure and thus the simple model used by Sliv and Band should agree 
well with real conditions. Green and Rose (1958) and Rose (1958) have re- 
cently recalculated the internal conversion coefficients taking the finite size of 
the nucleus into account, using a nuclear model differing from that of Sliv and 
Band. The so-called ‘‘no-penetration’” model of Rose confines the radiation 
sources to the nuclear center whereas Sliv and Band put them on the nuclear 
surface. In the no-penetration model, as in any physically reasonable model, 
the singularities at the origin are removed and thus we expect a general agree- 
ment between the results obtained from the surface model of Sliv and Band 
and from the no-penetration model as well as from any other physically reasonable 
model. This statement is justified by the investigations reported in the present 
paper. As a rule then, the contribution to the interaction between shell electrons 
and nuclear matter which comes from the intra-nuclear region is negligible as com- 
pared with the contribution from the extra-nuclear region. Provided that this 
is the case, the statement that the internal conversion coefficients are indepen- 
dent of the nuclear structure is still valid within a few per cent. The contribu- 
tion to the conversion matrix element from the extra-nuclear region is, however, 
proportional to the gamma-ray field and is thus very small when the gamma- 
ray emission is very weak. This contribution may then be of the same order 
of magnitude as the intra-nuclear contribution. In such cases, the internal con- 
version coefficients may depend considerably on the nuclear structure. Kramer 
(1956, 1957) has given the expressions for the conversion coefficients for a finite 
nucleus and Church and Weneser (1956) pointed out the possibility of large 
deviations of the conversion coefficients from the Sliv values and investigated 
some slow magnetic dipole transitions. Similar investigations have also been per- 
formed by Reiner (1958) and Nilsson and Rasmussen (1958). The latter authors 
investigated electric dipole transitions. 

In spite of the comprehensive tables for the conversion coefficients that are 
available, a theoretical investigation should be interesting both from the theo- 
retical and the practical points of view since such an investigation should, among 
other things, simplify the application of corrections for disturbing effects like the 
influence of the finite nuclear size. In the present paper the conversion process for 
energies near the conversion threshold is studied and the results are applied to mag- 
netic dipole conversion in the K-shell. The results obtained for the very low 
energies near the threshold hold also for energies which are as high as ten times 
the binding energy of the bound electron, and, for the heavy nuclei, extend 
well into the energy region of physical interest. 
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Since the detailed calculations leading to the expression for the internal con- 
version coefficients in terms of radial integrals may be found elsewhere, we 
will not attempt here a derivation of these expressions, but will adopt the results 
obtained by Rose et al. (1951) and Rose (1955). These results do not take 
the finite nuclear size into account. Such effects will be considered in the fol- 
lowing. 

In this paper, we shall confine ourselves mainly to conversion in the K-shell 
for magnetic transitions. For such conversion coefficients, the following expres- 
sion has been given by Rose 


pe = 2aak 
OR ee 


[L| Rvan+1(L) |? + (L+1)|Ryv=-1(L)/71, (1) 


where JL is the multipole order of the gamma-ray, « the fine-structure constant, 
k is the gamma-ray energy in units of mc”? and where R,,(Z) contains the radial 
integrals of the transition matrix elements. 

The expressions for the radial integrals become simpler and are more easily 
studied if a representation of the electron wave functions, introduced by one 
of us (Olsson 1959 a), is used. These expressions are 


Eut+y ¢ 
1= Za , H,t)+ @ +1,2,0)] 
Re Gains nerdy Gre ws 
Eut+y ¢ 
Xp = (x — lo , #, t) + ——__——_-® +18, 0], 
a= (%— +) (y ) Pea), (y ) 
where (2y),=2y(2y +1). 


Furthermore, Z is the nuclear charge number, H is the electron energy, in- 
cluding its rest mass, in units of mc’, t=2«Z r with r measured in units of 
h/me and finally y= +Vx®—0? Z? where x= +1, +2, +3, .... Both signs for y 
are possible and the functions thus obtained are linearly independent solutions 
of the radial part of the Dirac equation: 


df (r) aay yr (« 1422) gir 


dr r 
20) (p4.1+%2) f)-*** 909, y? 
dr * r 
where M%=rfl(r), X2=17g9(r). 
The functions O(y, #,t) are defined by 
D(y, H, t)=e"'" Fy (y tin; 2y; 2ipr) (4) 
where p=VE?-1 and n= as 
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For a point-nucleus, we must choose the positive value of y, otherwise the wave 
functions do not give finite electron densities at the origin. Introducing the 
wave functions (2) into the expressions for R,,(Z) as they are given by Rose, 
we obtain 

oe] 


Ry (L)= N,N; | hy (ker) (XP WP + UP AP] dr (5) 


0 


where the letters 7 and f refer to initial (bound) states and final (continuum) 
states, respectively. R depends on x’ of the final state. Angular selection rules 
forbid transitions other than those for which x’= —Z or +(£+1). The quantities 
N are normalization constants of the electron wave functions and h,(kr) is a 
spherical Hankel function of the first kind 


1 
1 (kr) = aes HY (k n=(-9' EF, (- La Aphids at (6) 
In agreement with common notations, the ,/y-function is defined by 


1 ign (eae BW Meg toy 
2Fo( - ad via 2, n! (se) 


where a symbol such as (a), means 


(2), =a(a+1)(a+2)... (at+tn—1), a=l. 


The expansion (7) breaks off when —Z+n—1=0 and is thus a polynomial 


of degree L in Dike 
at 

The spherical Hankel function can be expressed in terms of spherical Bessel 
functions. 


hy (ker) = jy (kr) — ¢(— 1) 7-1 (kr) (8) 


where jn (kr) = yest Tray (kr), 


Br IM hang pe J_ 1-4 (kr). 


This splits the Hankel function into its real and imaginary parts. Inserting 
(8) into (5) we obtain 


wo 


R,, (L) = N; Ni jn (kr) XOX? + AP WP] dr — 


0 
co) 


=i(= DENN, [ink OPH P+ UPaPidr. (9) 
0 
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The constants N; and N,, as well as both integrals, are all real. As is well-known, 
the Bessel function can be expressed in terms of confluent hypergeometric functions : 


1 er 
Ae 5 Rf re I - Dex), 
J, (x) para eh (vy +4; 2v+1; 2iz) 
We thus obtain 


ene 
in(kr)=3 (5) (E413 2L+2; 2ikn), 
2 


(10) 


aN a L+1 A. 
(-1) (=) D(E+) our pL, —2L; 2ikr). 


j-r(kr)= 2 rae Vin 


The last relation is not true if L is a non-negative integer but it can still be 
used provided that the integrals in which the function occurs in the integrand 
are defined also for the exceptional values of ZL. If the function itself is needed 


we can use 
; i (? +11 (L+4) Pe she 
j-r (kr) Wo Va of (2—-L; -—- 


but we prefer to use (10) here. 
The normalized wave functions for the K-shell are 


=e aZ(1—y), 2 
NiXy= Vere, é : 


ae aZ(1+y) y <5 
Nika | e202 : 


(11) 


were y=V1—o? 2. 
We can then express R,,(Z) in the following manner: 


Ry (L)=N; Vera) [ty ty a) Ja |x’ D+ 


"¢' +y’ 1 ’ ’ , , 12 272 
ee hy oy — A LL, 12 |. =Vx'*—o2Z. (12) 
oe @,'), y-Y a(Z, |x’ |), y 


Here all primed quantities refer to the final state that is uniquely deter- 
mined by x’ and its energy H’. The integral J, (L, |x'|) is 


Pai 
J, (Z, le = feo Fe ra Gen e Fi (y' +in; 2y’; 2iprydr (13) 
d 
: 
and J,(L,|x’|) is obtained by increasing y’ in J,(L,|x’|) by one unit. 


| 
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Inserting (10) into (13) we obtain for the real and imaginary parts: 
Re J, (L, |’ |)= 
arty’ IL V7 } ; 
= opty’ +L e ik’s 5 ie x 
Qa P (L+H) L+3) eis WEE 26425 20k 8) 


0 


oo 


eo 'PS Fy (y' +in; 2y'; 2ip’s)ds. (14) 
Im J, (L, |x’ |)= 


= aA TLAD) 0 ono gyty-t-1-ms w (_E; 2; 2488): 
ptt} Vx < : 
0 


-e 3 FL (y’ +in; 2y’; 2ip’s)ds. 


t l, 
Here, we have introduced 8=5, P=25 and v=. 
The two integrals (14) may be expressed in terms of the Appell-function 
F, (a, b, b’; c, c’; x,y) (Erdélyi 1936). 


grt KL Vx T(yt+y'+L+1) 
Re J, (L, |x’|)= ~ oy aE ; : 
DP Eep Deer (15) 
Qik’ Qin’ 
EF. ‘+E ik db a : ;2L 2, 2 es an -\? -770 ’ * 
aly+y Meares Hak nat leas Oe re ES aay 
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L, |x! 2 etre Ke 
Im J, ( > | x !) pert Vn [l+i(k’+p')y*”-* 
; Qik’ Qin’ 
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In the neighbourhood of the origin, the F,-function is defined by 


0 


ait (2) msn (b)n (Bn 2" Y" 
Tih DD 7 De OAc) eke, arate «RAN 
n=0 


Convergent expansions for integrals of the type occurring here have been 
given by Olsson (1958b). The expression for the imaginary part fails for 
L=0, 1, 2,3, .... We therefore introduce a transformation given by Erdélyi (1948). 


a 


i : 
F, (a, b, b'; 2b, ¢'3 2, y)= (1-3) Hi, (a, b+he's ro ie 54) (16) 


(A)em+n(B)n vy” 
=o (€)m(d)nm!n! ~ 


7Ms8 


where Ha, 05°C, a: ae 


re 
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Due to the introduction of the H,-function the J-integrals will be defined also 
for integral L-values: 
vty ttt Ve Ty t+y'+L4+1) 
L+1 A Vista ye J 
Re k J, (L, | x |) eae)": T(L+ 3) (ltip’)rty teri 


Hayy tO lay panel t8, 2y's2; y),20(12) 
Im kJ, (L, |x'|) = _ 2 (2a Z)'T(L +4) T(y+y'—L) 


Vx (14 ep yry-e 
“A, (yt+y'—L, y'+in; 4—L, 2y'; x,y), 
ry ide bs Qip' 
where x Bressit y ears 


By confluence in (16) we obtain 


1 
lim F(ab, —;2b,c’; vey) =F (a, 6; 2b, c'; 2, y)= 
e>0 é€ 
-—a 1 2 
=lim (1-3) Hy(a, 25044, epee et) (18) 


4 Wi ee 2Qy 
ad Se Jue oe 
(1 ;) H, (a; 644.0 42m) =), 


, = (2) min (0)m2™ y” 
h , : = 
where ye (a,b; ¢,¢; x, y) 2, (C)m (C")n m!\n! 


‘ ; = (B)emint™ Yy” 
and H, (a; 6, ¢; x, y) 2,, (b)m(C)n m!n! 
0 


These results are easily obtained by using the limit 


lim (« a; ‘) (ex)" = (ba)" 


e—>0 


in each term in the expansions of the functions in question. The parameter 
E2Z 
pipe tends to infinity when p—0O and the H,-functions in (17) will there- 


fore have its confluent function H, as a limit when p=0. Conservation of 
energy requires 
k=E-y 


and the gamma-ray energy for p=0 is then 
kj=1-y 
which is the threshold energy for conversion. 
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Internal conversion matrix elements near the conversion threshold 


From the result (18) of the previous section the threshold values of the 
J-integrals can be obtained directly: 


l-y 2L+1 Vn 
3) T (L+38) 
“Ho (yp +y' + L413 D+ b; By’ sas; Ue): 


Re kj** J, (L, | x’ 


)= 2 (2a.2)"( T(yt+y'+L+1)- 


1 
Im kJ, (L, [x)= - 27°" 2a y= nee UD yty'-L): 
3 


-Hy(yt+y'—L;4-L,2y'; x,y), (19) 


Ly. 


where Ly= — 4(1+y)’ 


Yo= —2. 


All transitions, the electric as well as the magnetic, can be expressed in 
terms of the J (LZ, |x'|)-integrals for which a convergent expansion can be ob- 
tained (Olsson 1959 c). This expansion is convergent for all values of the parame- 
ters, which are of interest in connection with internal conversion and has been used 
to some extent for checking results which have been obtained from simpler, 
approximate formulas. The complicated mathematical structure of the problem 
seems to make it impractical to try to obtain simple expressions for the con- 
version coefficients unless the gamma-ray energy is very low or very high. We 
will here focus our interest on the low energies. This simplifies the expressions 
considerably but the conclusions which can be drawn from the simpler for- 
mulas are in many cases of a very restricted validity since they are in general 
true only for very small energies. In the present work, we will investigate the 
internal conversion process in magnetic dipole transitions by taking advantage 
of the simplifications that are possible for low energies. 

Such an investigation is most interesting for the heavier nuclei for which also 
the validity of the simple expressions, which will be derived for the conversion 
coefficients, extend well into the energy region of physical interest. The energy 
dependence of the J-integrals is very complicated due to the fact that the 
energy is contained in the parameter 7 as well as in the arguments of the 
functions in the integrands and many terms seem to be needed in the conver- 
gent expansions of the integrals, even for very small p-values. However, if the 
validity of the expansion is restricted to very small p-values, an expansion in 
powers of p can be used and the result obtained using only a few terms may 
hold good far beyond the value of p where the convergence ceases. The pos- 
sibilities and short-comings of an analysis based solely on a low energy in- 
vestigation can be fairly well predicted from existing data. Besides, the in- 
formation obtained can be completed by an investigation of the high energy 
limit. 

In order to obtain an approximate formula for the conversion coefficients on 


j— 
the very threshold, we observe that 2) in (19) is small since Tat <0 for 
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all nuclei. The expansion given for the H,-function converges when |x|<4} for 
all values of y and the expansion of the ‘Ty, -function can thus be expected to 
be rapidly convergent when summed over x. The summation over y is easily 
performed. We put 


s (2)om (@+2mM)y,x™y" 


Oo EGS Til emer ar a 


=>) (@)2m x” . (a+2m)n -34 Phen 2" 


eS Gaml 2 eat Grnt Patt 2m: 03 9). 


For the small values of « obtained above, we simply put 
H, (a; b,c; «,y)~,F,(a;¢;y), |x|<} 
and obtain 


VOD (y ty + L+1) Va. 
T(L+® 


+1 we QrtY" ( 
Re ké*1J, (L, |x’|) 2a Z)! (52) 


iy ty POT; 2y 5 —2). 


C(yt+y' —D)T(L+3) 


Vx 


Im kj*? J, (LZ, |x’ |) ~ —27*” (2aZ) F,(y+y'-L; 2y'; -2). 


(20) 


In these expressions we have obviously neglected the binding energy of the 
converted electron. The most important contribution to the magnetic dipole 
conversion is the one where the converted electron is expelled in an s,-state. 
We then have x’= —1 and thus y’=y and hence we obtain 


2 e729” D (5 y(3—2y)(1—y)? 
Re KJ; (1, 1)» —3e°2”T (27) ee 2 
(21) 
Re KJ,(1, 1) © +2622” P (2y)2y (2y), 2 
ovg\> cms ee Leni 
It can be seen from (12) that the factor in front of the J, (1,1) integral in 


R_,(1) is (1+y)(2y), times larger than that in front of the J,(1, 1) integral. 
The real part of the latter integral is 2( (27) ),/(83-—2y) times larger than the 
real part of the former and is of opposite sign and hence there is a very good 
cancellation of the real terms. For the imaginary parts we obtain 


Im J, (1, 1) ~ —a Ze? 2" T (2y—1),F, (1; 2y; 2), 
Im kJ, (1,1) + —a Ze? 27471 (2y) F, (2, 2y +2, 2). (22) 


The terms which are omitted in (22) contain only ,F,-functions that are po- 
lynomials. In the expression for the J, (1,1) integral the error is of the order of 
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(l—y)? since the coefficient of (1—y) is zero. For light nuclei, the value of the 
imaginary part of the J,(1,1) integral is about half the value of the imaginary 
part of the J,(1,1) integral whereas, for the heaviest nuclei, the ratio is 
about 2. The two integrals are equal for a Z value of about 84. Furthermore 
we easily obtain 


Re J, (1, 1) 
Imy/5@5 1) 


Re J, (1, 1) 
Im J, (1, 1) 


| < 0.005, 


E 0.04. 


For light nuclei, these ratios are much smaller. Due to the above-mentioned 
cancellation of the real terms, the real part of the matrix element is much 
smaller than the above ratios would indicate. In fact, we have 


Re R_;(1) 
Im Ry (1) 


< 0.001. 


We can therefore safely neglect the real part of the matrix element. 
We now consider transitions to the ds),-states. A rough numerical estimate 
gives 


for all Z-values. Thus the contributions from the transitions to the ds;,-states 
can also be neglected for the very low energies considered here. 
With an accuracy of 10°-*, we can thus write at low energies: 


2 


VaZ(1—y) 
‘VP (+24) 


x Suak 


resi 3 


Jato A +e” ; 


== J g(a 23 
T= Beni] eB 


For the normalization constant N,, we take the expression (Olsson 1959 a) 


y’ , omy - ’ / 2 A 
x,=(5) E'xz IT (y' +in)|e (24) 
“Z x! — 2apl(1+2y) 


By means of asymptotic expansions of the [’-function, we obtain, for small 
p-values 


1 6x? y'+3xy'+y'—4y% 
r if be ue WA poo], (25) 


VaZ0 (1+2y’) 12 (x? — y") 


In order to obtain in addition the values of the J-integrals in the vicinity of 
the conversion threshold, we calculate their derivatives with respect to the 
energy neglecting the real parts. We obtain, using (14), 
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B+" (20 Z)ET (L+3]o)_ 
Vx 


et im J,(1, b= —- 


r I ajeere (yi tin; 2y'; 2tp’s) i of (3-2; Fre: 
0 


eal Be \ so _,, 
+ [eter etp, (4-2; = aa) ge Aly tins 24"; Zip's) ds}. 
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We first calculate the derivative of the Bessel function at the threshold 


r) k? s? s* b= 
= :' ube o ta ) > 2 5) 
lene (1 Hee reals erate eal _ et): a 


In order to obtain the derivative of the Coulomb wave functions, we use an 
expansion of these functions in terms of Bessel functions (Tricomi 1949, Olsson 
1959 a) 


e7'?* Fi (y' +in; 2y'; 2ip’s)= PA pe 
n=0(2')n 


(28)" Ff, (2y’ +n; —2Hs) (28) 


where the coefficients can be obtained from the recursion relation 


2 


(n+ laxza= — pep l(n ty’ — 1) dy-1— Bays], 


aj=1, a,=0 ete. 


Since the coefficient a, is of order of p*, only the four first terms have to be 
considered when carrying out the differentiation. The result is 


Dib) a>. aiiet es . os 
Fx i's Fy (y' +in; 2y'; 2p | 
mae i n+ 1, 2 helen. nal F,(2y'+2; —28)+ (29) 
tig © 1(2y ’ 8) a Z? (2y'+1)° wey ’ 
458° 
(Dy +8: — 28), 
+322 (2y'),° 1( yt 8) 


The derivatives calculated are inserted into (26). We observe that the integrals 
occurring will be of the same type as those which could be expressed in 
terms of H,-functions. After suitable changes of the parameters we obtain 
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_ PY Qa ZT (L+4) 


[Piyty'—L+2) 

(22-1) (1+) 
PS S| 

ee tO E, (pty —L+1)4 Dey bse, %)— (30) 
LGyty —L£+2) F a) Beotve a 

= 2D Qy +1) Ha (yy — L425 4D 2y' +25 x5, Ye) + 

40 (y+y'-L+3) 
3a? Z? (2y'), 


E Im k**? J, (1, | 


Hy (yt+y' —L+2;§—LD, 2y'; 2%, yo) — 


H(yry —L+334—L,.2y taste, w)| 


where X= 


as in (19). An approximate expression for the derivatives can be obtained by 
replacing the H,-functions by the first term in their expansions as was done 
in (20) and using recursion relations for the ,F,-functions to simplify the ex- 
pression obtained. 

With the aid of the results derived above, it is a matter of straightforward 
calculation to derive the coefficients in an expansion of the form 


36K = A,+(B-1)Bt+-. (31) 
We obtain Ag= (kB) z-1, 


B= 2A; E In eR.) 


E=1 
where 


a a 
In |? R_ sig eet 
E nf R (I) [ginw,| + 


op 4 E 2 ek he cedeeis 
(1+y) Ee Tat Ue ao PAWLV! +g ap, tral D 
. 


2 2 
(l+y) kod, Ch 2} 4 Qy), “0721s 1) 


+ 


For the derivatives we take the imaginary part (30). 
The expansion (31) is convergent for small values of #—1 only. Nevertheless 
the simple expression 


A+kB 
pi=-—5 (32) 


gives correct values of the magnetic dipole conversion coefficients for the 
heavier nuclei even for gamma-energies as high as 0,7 MeV. Since it is common to 
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Table 1. The coefficients A and B. 


Z | A | B | Z | A | B | thi 
30 0.000326 0.000244 64 0.0126 0.00927 4.10 
31 0.000375 0.000281 65 0.0138 0.0102 4.19 
32 0.000430 0.000323 66 0.0151 0.0111 4.27 
33 0.000492 0.000369 67 0.0166 0.0122 4.36 
34 0.000561 0.000420 68 0.0181 0.0133 4.46 
35 0.000638 0.000477 69 0.0199 0.0146 4.55 
36 0.000723 0.000541 70 0.0218 0.0160 4.65 
37 0.000817 0.000611 71 0.0239 0.0175 4.75 
38 0.000921 0.000688 72 0.0262 0.0192 4.85 
39 0.00104 0.000774 73 0.0287 0.0211 4.96 
40 0.00116 0.000868 74 0.0315 0.0231 5.07 
41 0.00130 0.000971 75 0.0345 0.0254 5.18 
42 0.00146 0.00109 76 0.0379 0.0278 5.30 
43 0.00162 0.00121 77 0.0417 0.0306 5.42 
44 0.00181 0.00135 78 0.0458 0.0336 5.55 
45 0.00201 0.00150 79 0.0504 0.0370 5.68 
46 0.00224 0.00166 80 0.0555 0.0407 5.82 
47 0.00248 0.00184 = | 81 0.0612 0.0449 5.97 
48 0.00275 0.00204 82 0.0675 0.0497 6.12 
49 0.00304 0.00226 83 0.0746 0.0547 6.28 
50 0.00336 0.00249 84 0.0825 0.0606 6.44 
51 0.00371 0.00275 85 0.0914 0.0671 6.62 
52 0.00409 0.00303 | 86 0.101 0.0744 6.80 
53 0.00451 0.00334 | 87 0.113 0.0826 7.00 
54 0.00496 0.00367 88 0.125 0.0919 7.20 
55 0.00546 0.00404 89 0.140 0.102 7.42 
56 0.00600 0.00444 90 0.156 0.114 7.66 
57 0.00659 0.00487 | 91 0.175 0.128 7.91 
58 0.00724 0.00535 | 92 0.197 0.144 8.18 
59 0.00795 0.00586 93 0.221 0.162 8.47 
60 0.00872 0.00643 94 0.250 0.183 8.78 
61 0.00956 0.00705 ) 95 0.283 0.207 9.12 
62 0.0105 0.00772 96 0.322 0.235 9.48 
63 0.0115 0.00846 97 0.367 0.267 9.88 
| 98 0.421 0.306 10.32 
99 0.485 0.352 10.81 
100 0.562 0.407 11.34 


refer to the gamma-energy, which is also the nuclear transition energy, rather 
than to the energy of the conversion electrons, we will use (32) instead of (31). 
The two expressions are simply connected by means of the relation k=H—y. 

The coefficients A and B have been computed for the range Z=30-— 100 and 
can be found in Table 1. By means of this table a very accurate and simple 
determination of the conversion coefficients ff can be made for any energy within 
the range of validity of (32). In order to give the reader an idea of the ap- 
plicability of (32), conversion coefficients derived from it and Table 1 are com- 
pared in Fig. 1 with values of k°f{ obtained by Rose. The values obtained by 
Rose have been complemented with a few exact calculations carried out by the 
present authors using the convergent expansions derived by one of them (Ols- 
son 1959 b). For the heavy nuclei, the agreement is good, even for energy values 
far beyond the domain of convergence, but for such values it is considered 
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Fig. 1b. 


Fig. 1. The magnetic dipole conversion coefficients for the K-shell, multiplied by k* and plotted 
as functions of k, for four different Z-values. The dotted lines are obtained from formula (32) 
in the text and should coincide with the point nucleus values computed by Rose. The two lower 
curves in each plot give the finite nuclear size values obtained by Sliv and Band and by Rose. 
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merely accidental and no attempt has been made to compute further coefficients 
n (31). 

Formula (32) gives the unscreened point-nucleus conversion coefficients but, 
since we are now convinced of its practical applicability, we can correct the 
coefficients A and B for various effects having an influence on the internal 
conversion process. The fact that both the constants refer to threshold energies, 
for which considerable simplifications in the mathematical formalism are possible, 
should make the investigation of such corrections as simple as possible. We 
will next consider corrections due to the finite nuclear size. Correction terms 
will be added to the electron wave function and, from them, the finite nuclear 
size wave functions in the vicinity of the nucleus can be obtained. 


Effects due to the finite size of the nucleus 


The wave functions given in (2) are valid without modification for the bound 
states as well as the continuum states. Furthermore two independent solutions 


are obtained from (2) simply by choosing for y either y= +Vx?—«?Z? or 
y= —V2—0?Z. In order to obtain values for the wave functions in the vicinity 
of the origin, we can use the power series expansion (Olsson 1959 a) 


Oy, H,)= > 0 (33) 


Cn 
gs 


where the coefficient can be derived from the recurrence relation 
P 


(1 +24) Cn41+ Bey + 22 


| _2 72 Uen- — =( 
Co=1, c-n=0 


We shall first derive the wave functions for bound states. Since the presence 
of the nucleus no longer makes it necessary to require %,(0)=%,(0)=0, the 
wave functions will contain both the regular and irregular solutions of (3) but 
they are still submitted to the condition of normalizability. Normalization is 
possible only if the solutions decrease exponentially at infinity. The correct 
mixing of the regular and irregular solutions can be obtained from the well-known 
asymptotic behaviour of the wave functions and the result, which contains four 
® (y, H,t) functions, can be conveniently expressed in terms of Whittaker func- 
tions. 


n H(x—y) 
t= 0 | Wasn a Oy Bigay Narnren a 
n E(x+y) 
hy =(%—y) Warnes SE RtOE Ea Tea Wriyytt o) 
where e=2V1-H’r, n= LEE -». 
Vli-E? 
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The two Whittaker functions occurring here are defined as certain linear com- 
binations of confluent hypergeometric functions: 


| Res ee ea: 2y; ")+ 
ao girth e #2 FB (l—n—-2y; 2—2y; 2) 
35 
Wasyyey aay att et FB (l—n; 2y +2; x) + oe 
+ Tan event pe (—a— 25 — 2; 2), 


If the parameter n is a positive integer or zero, the irregular terms in (35) 
are zero and both the components 7%, and Z, vanish for r=0. This gives the 
ordinary hydrogen-like wave functions. In the presence of a nucleus of finite 
dimensions, the exterior solution (34) must be joined continuously to an interior 
solution at some point outside the nuclear charge distribution. This can be ac- 
complished by changing the parameter n which is then no longer an integer. 
Due to the small dimensions of the nucleus, the deviation of n from an integral 
value is very small and sufficiently accurate wave functions can be obtained 
by considering only terms of first order in the deviations 6. Since solutions 
that are regular at the origin can be obtained for integral n-values only, it fol- 
lows that the terms of first order in 6m are singular at the origin and thus 
considerable changes in the electron densities can result in spite of the small 
value of 6m but only where the terms containing r” are large, i.e. for 
ram: 

Where r>1, the corrections are of order 6m and are thus negligible. Since 
we are here concerned with K-shell wave functions only, obtained from (34) if 
we put n=0, E=y and x=-—1, we shall investigate (34) when n differs by 
a small amount dn from zero. We obtain, retaining first order terms only, 


L+ 
L,=aZ | Wonsyya (2) Snag? Wines (| 
rd! (36) 
Ager (ich y) | Woneny-s (2) bmn Wiepts o) : 


From these expansions it is immediately seen that the small component, %,, is 
much more sensitive to the presence of nuclear charges than is the large com- 
ponent. We are only interested in the terms in W,,,,,(«) which can be large 
and thus we put 


Dla 2). « 
Wii’, (aywr OTA), 2 


This is exactly equal to the second term in W,,,,, (a), as can be seen from 
(35) since ,F,(—2y; —2y; x)=e". Similarly we get 
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K-shell wavefunctions for Z=82 
in the vicinity of the nucleus 


Fitted to constant charge 
distribution 


0 
}-X,  én= 0.001 
0.002 


3 -13 
nuclear radius (1.3°¥A -10 cm) 


Fig. 2. Changes in the small and large components in the electron wave function near the 
nucleus due to variation of On. The functions are obtained from formula (37) in the text. 


The small component is plotted, for convenience, with an opposite sign. 
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For the sake of simplicity we write 


_2p_6nT (2y-1) 
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which finally gives the K-shell wave functions 
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(37) 


The corrected wave functions for reasonable values of dn are given in fig. 2. 

It should be mentioned that the change in the normalization constant is of 
order 6” and can thus be neglected. The value of 6 is obtained by joining the 
solutions (37) to interior solutions at some properly chosen point x) which may 
be taken as the nuclear radius. The value of the quotient of the interior solutions, 
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(; a) interior 
Xz (%q) : 


depends on the nuclear charge distribution but it is very insensitive to variations 
in the energy due to the small value of the nuclear radius. This quotient de- 
termines 6” which is obtained from 


én ” Xn 


Xz (Xo) ty 3 on Ss lens 
1 D (27) 23 (l—y)e aes 


(38) 


interior 
Since (2) is negative for the small values of 2) considered here, we 
2. 


must have 


én i La 
1 Dy) s Jaen Pras >0 


xe” 


D(2y) la cers 4 


or on< 


which gives 6n<0.0016 even for the heaviest nuclei. 

The continuum wave functions will now be treated in a similar manner. 
We shall use %* to denote the solutions (2) when y is positive and Z% the 
solutions when y is negative. The normalization constant for the X%~ wave 
functions is obtained from (24) by changing the sign of y. 

The normalized free wave tunctions with an arbitrary asymptotic phase can 
be written 


Nt one 
ag is Me el berg 

N} Nz em) 
Ser + bm, te 


where A and 6n,; are two parameters. A depends on 6m, since only one para- 
meter remains in the general solution if the solutions are normalized. The nor- 
malization constants Nj}, chosen as in (24), give the following asymptotic ex- 
pressions 


Xt © / 


— sin prin In 2prz —” —arg P (by + in) + "| 
ap 2 


(40) 
We 0 Ae ae |pren In 2p rt — arg P( by + in) +3" 
X 
Z(H-1 
where tan 6* a ca a 
p(x+y) 
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Generally we have 


[isin (w+ 6,) + 6m, sin (w+ 6,)]=sin (2+ 63) 


- [cos (w+ 6,) +6, cos («+ 6,)] = cos («+ 63) 


_ sin 6,+6n,; sin 6, 


if tan 6 
° cos 6, +6 n, Cos dy 


and A? =1+426n,; cos (6, — 65) + (6n,)*. 


The solution (39) thus has the asymptotic form (40) but whith a different 
phase. Further it is seen that the changes in the normalization constants are of 
the order dn, which, for small kinetic energies of the free electron, is of the 
same order as 6n, introduced into the wave functions for the bound states. 

Therefore, only a few terms in the expansions of % are needed and these 
terms can be obtained from (33). We obtain 


aZ Lit. "Bee. ot 
1 Ae 
* tae “op, es | 


a ae Et Ex-y t 
a 1 - 
tae ES OMe? 


and, for the free electrons, we thus obtain 


OL, 1p eee 1 
=N} Xt +6n,N7 1 ( + ) -.| 
MEA LAE ov | Ney ey (—2ye) 


“@+y) E Ex-y 1 
= Nt yt aA SA 
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The parameter dn, is obtained by connecting the solutions (41) with the 
interior solutions as in (38). 
For an s,-state, for which x= —1, we have 


Z 
y= Nj Ht + 9mpN7 AF [1 : (z+, Et+y Jen] 


(41) 


2 l+y)(2y-1 
eas y ls ) (42) 
mpeg Ati Lay t thee ak oe” bas 
X= Nz Xz — On, Nz t” Re ee | 


The correction terms should, of course, be used only for t¢<1. 

The most important corrections due to the finite nuclear size come, how- 
ever, from the effects inside the nucleus where the rather strong singularities 
of the fields (which nearly make the integrals non-existing (1)) considerably ex- 


(1) In the limit of a point nucleus the integrals may actually diverge if the limit is not 
carefully derived. This circumstance has given rise to the so-called gauge-problem, the spurious 
nature of which was first pointed out by Kramer (1956, 1957). 
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aggerate the interaction between the shell electrons and the nucleus. If the 
nucleus is extended to finite dimensions, both the photon-field and the electron 
wave functions are finite everywhere and the interaction between the nucleons 
and the electrons is more correctly estimated when the latter are inside the 
nucleus. To neglect the interaction completely in this region (Sliv) is, physically, 
a more reasonable approximation than to consider the nucleus as a point. 
However, in transitions where the integrand of the matrix element is not 
singular, or only weakly singular, the point-nucleus approximation is good. 
The latter statement holds for the majority of nuclear transitions. An im- 
portant improvement of the theory of internal conversion was thus made by 
Sliv who removed the incorrectly retained singularities from the theory. That 
this removal of the singularities is the essential effect will become clear from 
the present investigation. Indications in the same direction are also evident 
from the calculations of Rose (1958). This author used a model where the radia- 
tion sources were concentrated to the origin (no-penetration model). As a result 
he had the strong near-zone radiation field inside and in the vicinity of the 
nucleus. In this way, the strength of the interaction near the nucleus will be 
exaggerated but, nevertheless, only small differences from the numerical results 
of Sliv and Band were obtained. The reason for this is the use of finite nuc- 
lear size electron wave functions which removes the singularities in the matrix 
elements. Thus, for gamma-transitions that proceed at a normal rate, the con- 
version coefficients are insensitive to the nuclear structure. The differences bet- 
wen the calculations of Sliv and Band and those of Rose indicate how large a 
dependence on the nuclear structure that can be expected. 

A detailed analysis of the effects referred to above can be made by investi- 
gating their influence on the parameters A and B in formula (32). For the 
heavy nuclei, for which the applicability of (32) extends well into the energy 
region of physical interest, sufficiently accurate results can be obtained for any 
nuclear model.1 We will not investigate here the influence on the internal con- 
version coefficients that can be expected from the nuclear structure effects. In 
order to do this for a definite nuclear transition a particular nuclear model and 
the nuclear wave functions for that particular model and nuclear transition are 
required. Such investigations are only motivated for very slow gamma-transitions 
which are many orders of magnitude slower than the single particle estimates. 
Not knowing in detail what to put into the nucleus we shall thus use an empty 
nucleus in the same way as Sliv and Band. When necessary, the empty nucleus 
can be replaced by suitable expressions which take account of the nuclear struc- 
ture and the conversion coefficients can be corrected for the effects of nuclear 
structure in a way similar to that of Church and Weneser (1956). 

We will therefore not integrate over the nucleus and the corrected matrix 
element can be written 


oe 


Ry (L)= —1(—1) NN, fi-. (kr) (XA te + x1 x2) dr, (43) 


To 


where the electron wave functions are now given in (36) and (41). In order to 
obtain an expression simple enough for the correction term in R,,(Z) we retain 
only the leading term in the corrections added to the wave functions. 


+ Strong deviations from the spherical shape may form exceptions. 


381 


P. O. M. OLSSON, S, HULTBERG, A theoretical study of the internal conversion 


2 pe _ 72 are |/ eZ (1) hes 
2 Re, (1) = By (1) +i. Nj 220-0) 2(L+y) | S25") art 
0 


ee Ye 1 
x, aZ(l—y) f dr 

+210N; ater a+r eyon—-Ze dn, {g- 

oe 
P(1+2y) vy Phase | 


=k* R_,(1)+iN} sain 2 


al 
+ [Penden 2m 75 | *}. (44) 


The first term is the point-nucleus matrix element and the second term con- 
tains the corrections due to the finite nuclear size. The term containing 19’? 
is the most important one and removes the unduly strong interaction connected 
with the point-nucleus model. The terms containing the parameters én and 6 ny, 
represent the static corrections which arise from effects due to the changed electron 
densities near the nucleus. The nuclear radius, 7), is measured in units of h/mce. 
The two parameters 6” and 6m, are not independent but connected through 


the approximate relation 

x final state x initial state 

Gs et hee wo 
Xo Te Xo T% 


since the two ratios have to be equal to the ratio of the interior solutions 
which is insensitive to changes in the energy. From (45) it follows that 


ha  bnp~ —(1+y) TD (2y) dn. (46) 
Nj 


Writing (44) as IR, (1)=h Bi (1) +A R_1 (1) (47) 
and using (46) we obtain 


aZ(1—y) 


1 21 (2) dn 
T(1+2y) 
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The correction term is energy-dependent only in the factor N, which here is 


1 
Ne Y 2 
ony Faamonld ay ae . 

Then E RAR, (| = =7 [2 A R-4(1)]e-1- (49) 


We now write the internal conversion coefficients, corrected for the finite 
nuclear size, in the same way as that used in (31) 
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= "2% (1x R|—|k2 A RIP = 


= A+ (H-1) Bo++--=A°C+ kh B+. 
where A°= AZ— B°(1—y). (50) 


For simplicity, we have dropped the indices for R. We obtain 


ARI\? 
A Rl 
and finally B= (1 ae ) E <AsY As IF (52) 


The corrections to be made in A and B thus can be expressed in terms of 
the ratio AR/R for which the general expression is 


AR| 4aZ(1+y) (2aZr,)"} 2 T(2y)dn 
Pe a rey ay P42@r— edn ie 
where 7’ is given in table 1. The general expression is 
2 
T=2(1+y) kJ, (1, 1)+ 2») kode lll): (54) 
rt 2 


In (53) we observe that the second term contains the factor I’ (2 y)dn/(2aZ 19)" 
which is a quantity that can be determined from the eigenvalue equation (38). 
The approximate value of 5 was discussed in connection with (38) and com- 
paring with (53) we can easily see that the second term in (53) is rather small 
compared with the first term. 


Magnetic dipole conversion coefficients with finite nuclear size corrections 


In this section, we shall apply the results of the previous sections to mag- 
netic dipole K-shell conversion in lead. From Table 1 we obtain, for the point- 
nucleus, the conversion coefficient as function of the nuclear transition energy k: 


1 
pr= B (0.0675 + 0.0497 k). (55) 
_ This simple formula gives the ff obtained by Rose with an error less than 
me per cent for transition energies up to 600 keV. Further we have 
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Fig. 3. The ratio at the nuclear surface between the small and the large components of the 
interior electron wave function. This ratio is very insensitive to the choice of the nuclear 
radius and to the changes in the energy. The ratios plotted here are computed for r,=0 and H=1. 


| A A _ 2.894 x 1.801 (1.197 79)° 


Z saan aI [1 +2x 0.602 a] (56) 


where a2 can be obtained from the charge distribution within the lead nucleus, 
since 

T(2y) dn 

(2a Zr)” 


The quantity x can be obtained from formula (38) and the value of %,/%, from 
Fig. 3. For Z=82, we obtain 2~0.15 and then the correction due to the change 
of the electron density outside the nucleus will be about 18 % of the correction 
due to the removal of the unduly strong fields inside the nucleus. Fig. 4 shows 
the conversion coefficient for two different values of the nuclear radius, for a 
constant charge distribution within the nucleus, and for the nuclear charge 
distributed over the nuclear surface. The results are in good agreement with 
the values obtained by Sliv and Band. The small differences are due to the 
use of the approximate formula (53) and the neglection of screening in the 
present work. Screening is thus, at least for small energies, unimportant for 
magnetic dipole conversion in the K-shell. 

Since both the coefficients A, and B, defined in (31), are insensitive to the 
value of the binding energy of the converted electron, the coefficient A, defined 
by A=A;—B(1—y), can easily be corrected for the true binding energy H;,. 
We then have 

A=A,;—B Ez,. 


This increases the values of the conversion coefficients for the heavy nuclei 
by one or two per cent. For lead the values of Fig. 4 should be increased by about 
one per cent. The error in (53) is, however, perhaps two per cent since only the 
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Fig. 4. Magnetic dipole conversion coefficients for conversion in the K-shell, computed for 
different nuclear charge distributions and for different values of nuclear radius. The dotted 
lines represent results obtained by the present authors. The curve valid for the smaller 
nuclear radius and for constant charge distribution should agree well with the values ob- 
tained by Sliv and Band since the essential difference between the two calculations is that the 
latter authors have taken screening into account. Screening should not be of importance in 
this case since the transition matrix elements are only sensitive to the electron wave func- 
tions near the nucleus. The two upper curves represent point-nucleus values. 


leading term in the expansions of the wave functions has been retained, in 
order to obtain a simple expression. From Fig. 4 it is clear that distinctions 
between possible values for the nuclear radius and reasonable nuclear charge distri- 
butions require an experimental accuracy of a few per cent. A few careful absolute 
measurements of M1 conversion coefficients are shown in Fig. 4 (Stockendal and 
Hultberg, 1959). It should perhaps be pointed out that these measurements do 
not represent the best accuracy possible with the method of Hultberg and Stoc- 
kendal (Hultberg and Stockendal, 1959). In favourable cases an accuracy of better 
than 5% may be reached. As far as these experimentally studied transitions can 
be regarded to have a pure M1 character the preference for the finite nuclear 
size values can be seen. To within a few per cent, the latter are obtained by 
removing the singularity at the origin which, as has been stressed previously, 
greatly exaggerates the interaction between shell electrons and nucleons when 
the electrons are in the vicinity of the nucleus. The remaining corrections that 
depend on screening, the nuclear charge distribution, deviations from the Coulomb 
potential near the nucleus due to non-spherical shape of the nucleus and the 
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radiation mechanism inside the nucleus can amount to about 12 per cent as 
can be seen from Fig. 4. For nuclei, heavier than lead, these corrections are 
larger. Since the static finite nuclear size corrections depend more on the nu- 
clear radius than on the variations of the nuclear charge distribution, the models 
investigated here are not unreasonable. It is therefore concluded that, even for 
unhindered nuclear transitions, deviations of the order of 10 per cent from the 
values calculated by Sliv and Band cannot be excluded. Such deviations should 
be particularly interesting for deformed nuclei. 
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